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Abstract 

Three  model  reduction  techniques  are  examined  for 
their  applicability  to  large  lightly  damped  systems. 
Frequency  truncation,  modal  cost  analysis,  and  approxi¬ 
mate  internal  balancing  methods  are  applied  to  the  CSDL  I 
and  II  models.  Two  actuator/sensor  configurations  are 
investigated  on  the  CSDL  I.  Control  of  the  CSDL  I  model 
is  implemented  through  the  use  of  steady  state  optimal 
regulator  theory,  and  the  effect  of  various  reduced 
order  control  models  on  the  structure's  time  response  is 
detailed.  Modal  cost  analysis  and  internal  balancing 
yield  essentially  equivalent  results  for  the  cases 
chosen,  although  the  internal  balancing  technique  pro¬ 
vides  more  information  concerning  actuator/sensor 
effectiveness.  These  methods  provide  more  effective 
reduced  order  models  for  higher  order  systems  than 
simple  frequency  truncation. 


MODEL  REDUCTION  TECHNIQUES  APPLIED  TO 
THE  CONTROL  OF  LARGE  SPACE  STRUCTURES 


I.  INTRODUCTION 


The  current  state  of  the  art  in  space  flight  includes 
the  ability  to  orbit  and  manipulate  large  cargos  via  the 
space  shuttle  system.  This  capability  makes  feasible  the 
construction  of  large  structures  while  in  orbit;  these 
structures  may  take  the  form  of  space  stations  (for  commer¬ 
cial,  scientific,  and  military  use) ,  large  communication 
antennas,  and  orbiting  observatories,  to  name  a  few  possi¬ 
bilities.  Two  common  characteristics  of  many  of  these 
structures  are: 

1. )  use  of  low  mass  structural  elements  -  weight  and 

size  restrictions  of  the  payloads  will  dictate 
structures  designed  with  thin,  lightweight  members, 

2. )  low  damping  -  as  a  consequence  of  the  above  con¬ 

straints,  the  structures  will  exhibit  low  damping 
to  vibrations,  especially  as  the  size  of  the  struc¬ 
ture  increases. 

Regardless  of  their  highly  flexible  nature,  these  large 
space  structures  will  be  required  to  support  activities  depend¬ 
ent  on  stable  platforms.  Antennas,  lasers,  and  large  optical 
telescopes  are  a  few  examples  of  systems  which  will  be  depend¬ 
ent  on  precise  pointing  ability  to  perform  satisfactorily.  In 
general,  passive  damping  achieved  by  judicious  material  selec¬ 
tion  and  construction  techniques  will  not  be  sufficient  to 
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achieve  this  stringent  structural  stability.  The  problem 
facing  the  control  designer,  then,  is  to  damp  the  unwanted 
vibrational  motion  while  meeting  the  required  performance 
objectives . 

To  achieve  the  above  results,  active  feedback  control 
methods  are  required.  The  structure  is  first  modeled  mathe¬ 
matically,  and  its  structural  modes  and  mode  shapes  are 
identified  from  the  solution  of  the  undamped,  free  vibrational 
problem: 

Mq  +  Kq  +  0  (1) 

Finite  element  methods  (NASTRAN)  are  utilized  in  this  study. 
Control  is  achieved  through  the  use  of  active  control  tech¬ 
niques  appled  to  a  reduced  order  model.  The  issue  of  imple¬ 
menting  the  required  force  actuators  is  not  addressed  here; 
they  are  assumed  to  exist  in  a  functional  form.  Measurements 
of  the  success  of  the  control  are  made  by  strategically  placed 
position  and  velocity  sensors.  Using  this  sensor  output,  a 
feedback  control  law  is  generated.  See  section  III  for  fur¬ 
ther  discussion. 

Model  Reduction  Problem 

The  above  control  problem  is  very  large.  Any  real 
structure  has  an  infinite  number  of  vibrational  modes,  of 
which  only  a  relative  few  are  modeled  mathematically.  Even 
after  this  inherent  truncation  of  unmodeled  modes,  the 
control  of  all  modeled  modes  is  generally  an  intractable 
problem;  hundreds  of  modes  may  still  remain.  The  computa- 
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tional  burden  of  controlling  such  a  system  would  be  very 
high,  and  its  on-line  use  would  be  beyond  the  capabilities 
of  even  the  present  state  of  the  art  computers.  In  addi¬ 
tion,  it  is  unnecessary  to  control  modes  not  significant 
to  the  performance  objectives. 

It  is  thus  desireable  to  have  some  scheme  to  reduce 
the  control  model  in  size.  However,  the  reduced  model  must 
retain  the  significant  dynamical  characteristics  of  the 
entire  structure,  and  give  a  true  picture  of  the  control 
effectiveness  in  achieving  the  desired  performance.  This 
thesis  considers  selective  deletion  of  modes  from  the 
control  problem  while  still  recognizing  their  impact  on 
the  full  mathematical  model  of  the  structure.  In  effect, 

”a  selected  set  of  the  differential  equations  of  the 
system  model  are  discarded"  (18) . 

If  the  input-output  relation  of  the  system  can  be 
expressed  in  terms  of  a  transfer  function,  the  deletion  of 
modes  corresponds  to  eliminating  poles.  Retained  poles  are 
characterized  as  dominant  (those  closest  to  the  imaginary 
axis) ,  and  the  transfer  function  of  the  reduced  order  model 
will  ideally  represent  that  of  the  full  model  with  as  small 
an  error  as  possible.  However,  the  representation  of  a 
large  scale  multi-input/multi-output  system  by  matrices  of 
transfer  functions  is  unwieldy;  state  space  methods  are 
typically  used,  and  will  be  employed  here. 


REAL  WORLD 

MATHEMATICAL  MODEL 

DESIGN  MODEL 

IMPLEMENTED  MODEL 

FIGURE  1.  Model  Reduction  Process 

Levels  of  Reduction 

The  model  reduction  problem  is  characterized 
in  Figure  1.  An  actual  structure  has  an  infinite  number 
of  vibrational  modes.  The  first  reduction  occurs  in  the 
process  of  modeling  the  structure  mathematically.  Finite 
element  techniques  represent  the  actual  structure  as  a 
collection  of  point  masses  and  basic  structural  elements, 
often  with  ideal  characteristics.  Here  the  infinite  degree 
of  freedom  system  is  reduced  to  a  finite  degree  of  freedom 
system.  Although  the  subject  of  some  controversy,  the 
assumption  in  this  paper  is  that  the  NASTRAN  model  is  a 
good  representation  of  the  true  system. 

A  further  step  is  the  reduction  of  this  mathematical 
model  to  a  tractable  size  for  computational  purposes.  In 
this  step  it  must  be  determined  which  of  the  modeled  modes 
are  most  significant  to  the  control  problem  at  hand.  The 
retained  modes  at  this  point  constitute  the  design  model. 
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After  the  design  model  control  problem  is  solved, 
the  solution  may  still  be  too  unwieldy  to  implement  in 
real  time  on  an  actual  system.  Further  truncation  of 
the  closed  loop  system  may  then  be  necessary. 

This  thesis  addresses  the  model  reduction  process 
from  the  finite  element  mathematical  model  to  the  design 
model.  The  intent  is  to  identify  the  modes  to  be  active¬ 
ly  controlled  in  order  that  the  performance  objectives 
are  met  by  the  full  system  mode .  Three  approaches  to  the 
model  reduction  problem  are  examined  and  applied  to  the 
structures  described  below.  Based  on  these  reduced  models, 
active  feedback  control  is  then  utilized  to  achieve  the 
required  pointing  objectives,  and  differences  due  to  the 
reduction  methods  are  explored. 

Description  of  Models  Used 

In  order  to  facilitate  comparisons  between  various 
controller  design  methods,  several  standardized  structures 
have  been  developed  by  the  Charles  Stark  Draper  Laboratory. 
Two  of  these  models  are  utilized  in  this  thesis. 

The  CSDL  I  model  (Figure  2)  is  a  basic  tetra¬ 
hedral  structure  supported  by  six  legs  attached  to  ground. 
Such  a  tetrahedral  arrangement  is  considered  to  be  a  funda¬ 
mental  building  block  for  larger  assemblies.  The  NASTRAN 
analysis  utilizes  ideal  massless  truss  members  which  trans¬ 
mit  axial  forces  only;  no  torsional  or  bending  moments  are 
modeled.  Point  masses  are  lumped  at  grid  points  one  through 


four,  the  intersections  of  the  tetrahedral  arms.  Each  mass 
exhibits  three  degrees  of  freedom,  yielding  a  total  of 
twelve  for  the  structure.  Six  co-located  actuator/sensor 
pairs  are  nominally  located  on  the  structure's  legs.  The 
performance  index  is  considered  to  be  the  displacement 
of  node  1,  which  simulates  a  pointing  line-of-sight  require¬ 
ment.  The  typical  criterion  for  a  successful  control  de¬ 
sign  is  to  reduce  the  x  and  y  displacements  to  less  than 
0.0004  radians  and  0.00025  radians  respectively  in  20  seconds, 
given  specified  initial  conditions  with  no  additional  dis¬ 
turbances  .  This  study  uses  the  displaced  radius  from  a 
zero  position  of  node  1  for  time  response  comparisons. 

The  CSDL  II  model  (Figure  3)  represents  a  non¬ 
trivial  approximation  to  a  three  mirror  optical  space 
system.  Structural  members  include  bending,  axial,  and 
torsional  stiffnesses  in  the  finite  element  model.  The 
total  structure  is  nearly  28  meters  high  and  has  a  mass 
of  9300  kg;  twenty-one  co-located  actuators  and  sensors 
and  two  sinusoidal  disturbances  complete  the  model.  This 
model  was  used  in  a  multiple  controller  investigation  by 
Calico  and  Aldridge  (1)  and  in  an  application  of  direct 
output  feedback  methods  by  Calico  and  Thyfault  (23) .  Here 
it  is  used  primarily  to  exercise  the  model  reduction  schemes 
of  interest  for  a  large  scale  system.  Control  implementa¬ 
tion  and  system  time  responses  were  not  achieved. 

Further  configuration  details  for  the  two  models  and  the 
major  NASTRAN  analysis  results  are  given  in  Appendix  A. 


Q  6/6  co-located 
act/sens  case 

4  3  actuator  case 
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II.  EQUATIONS  OF  MOTION 


The  following  relationships  hold  for  the  systems 
of  interest  in  this  study.  For  a  large  structure,  the 
vibrational  equations  of  motion  can  be  expressed  as 

Mq  +  Eq  +  Kq  =  Du  +  Rw  (2) 

where 

M  =  nxn  mass  matrix  R  =  rxr  disturbance  map 

matrix 

E  =  nxn  damping  matrix 

K  =  nxn  stiffness  matrix 

D  =  nxm  actuator  map  matrix 

q  s  nxl  physical  coordinate  vector 

u  5  mxl  control  force  vector 

w  =  rxl  disturbance  vector  with  statistics 

E(w(t)}  =  0  E{w(t)w(r)T}  =  W6 ( t-T )  (3) 

W  =  rxr  noise  strength  matrix 

Solving  the  undamped  free  vibration  problem 

Mq  +  Kq  *  0  (4) 

yields  the  system  natural  frequencies  ok,  i=l...n,  and  the 
nxn  modal  matrix  i  such  that 

q  -  4>n  (5) 

where  q  is  the  nxl  vector  of  modal  coordinates  and  the 
columns  of  $  are  the  eigenvectors  u  of  [m&»?  -  K  =  0. 


T 


Substituting  (5)  i-n  (2)  and  premultiplying  by  $  gives 


-  £tdu  +  £trw  (6) 


which,  due  to  the  properties  of  <f,  is 

-  "1  i— 

\  |\  T  T 

n  +  n  +  I  to?  n  =  0  Du  +  $  Rw 


L  N  L  \l 

Note  that  the  damped  case  does  not  in  general  present  a 
diagonal  damping  matrix  in  modal  coordinates.  However, 
when  the  damping  values  are  small,  as  they  are  here, 
the  off-diagonal  terms  can  be  considered  as  second-order 
terms  and  a  diagonal  approximation  may  be  achieved  by 
discarding  them  (12) . 

Now  consider  an  "output"  or  index  of  performance. 


y  =  C  q  +  C  d 
i  _p2  -v- 


where 


0^  =  pxn  displacement  matrix 
Cv  =  vxn  velocity  matrix 

such  that  y  is  some  linear  combination  of  states  of  interest, 
In  the  problem  at  hand,  y  represents  the  line  of  sight  de¬ 
viation  from  nominal. 


A  similar  relation  is  the  measurement  vector  of  sensor 


output 


l  -  ?p2  * 


where  H  and  H„  are  position  and  velocity  sensor  map  matrices 

-p  —V 


'•  V-V- 
.  *  '.1 


Again,  substitution  of  (5)  yields 


y 


c  4>  n  + 

-p-j  - 


n 


(10) 


2  =  Hp$  n  +  Hv$  n 


In  state  space  form,  the  system  equations  become 


(12) 

(13) 

(14) 


x  =  Ax  +  Bu  +  Gw 
y  =  Cx 
z  =  Hx 


H  = 


(17) 


These  equations  form  the  basis  for  the  following  development. 


III.  MODEL  REDUCTION  METHODS 


Many  model  reduction  methods  exist,  and  are  typically 
grouped  into  the  categories  of  parameter  optimization, 
polynomial  approximation,  and  component  truncation  (18). 

In  the  first,  a  quadratic  error  criterion  is  minimized  by 
determination  of  the  parameters  of  a  system  through  non¬ 
linear  mathematical  programming.  The  second  group  attempts 
to  approximate  the  system's  frequency  response  over  a  given 
range  of  frequencies ,  or  alternatively  tries  to  truncate 
an  expanded  version  of  the  matrix  transfer  function  to  arrive 
at  a  reduced  model.  The  third  category  endeavors  to  retain 
a  select  set  of  the  differential  equations  of  the  system 
(the  "components")  while  truncating  the  rest.  Three  methods 
that  belong  to  this  third  group  are  frequency  truncation, 
modal  cost  analysis,  and  internal  balancing;  these  tech¬ 
niques  will  be  investigated  here.  Their  selection  was  based 
on  1)  ease  of  applicability  to  large  scale  multiple  input, 
multiple  output  systems,  and  2)  compatibility  with  the  finite 
element  representation  of  the  structures  of  interest. 

Frequency  Truncation 

The  simplest  method  involves  merely  truncating  the 
higher  order  modes.  For  many  control  problems  this  is 
adequate;  the  low  frequency  modes  contribute  most  to  the 
performance  objectives,  and  the  higher  frequency  modes  are 
essentially  outside  the  required  response  bandwidth.  An 
example  is  an  aircraft  control  problem  that  may  require  a 


designer  to  account  for  effects  due  to  low  frequency 
bending  modes;  higher  modes,  through,  simply  do  not  impact 
the  control  problem. 

A  different  situation  exists  for  large  lightly  damped 
structures.  First,  these  assemblies  are  much  less  stiff 
as  a  whole;  finite  element  analysis  reveals  a  cluster 
of  vibrational  modes  within  the  bandwidth  of  the  control¬ 
lers  and  sensors.  Second,  the  pointing  requirements 
of  onboard  systems  may  be  much  more  stringent,  making 
the  consideration  of  the  effects  of  higher  modes  impera¬ 
tive.  A  simple  truncation  may  still  leave  a  large  number 
of  modes  to  be  dealt  with.  Even  so,  one  tends  to  favor 
the  retention  of  lower  frequency  modes  when  possible 
since  the  model  is  known  more  accurately  at  low  frequen¬ 
cies,  a  lower  sampling  rate  on  the  control  law  will  be 
possible,  and  the  actuators  and  sensors  do  have  a  finite 
bandwidth  (16) . 

Modal  Cost  Analysis 

Modal  cost  analysis  is  a  model  reduction  technique 
developed  principally  by  Skelton  (7,18,19,20,21)  and 
is  a  special  application  of  his  more  general  component 
cost  analysis  (18).  Essentially,  a  cost  function  is  es¬ 
tablished  detailing  the  weighted  contributions  of  each 
component  of  a  system  to  an  output  relationship.  The  indi¬ 
vidual  contributions  of  each  component  are  assessed,  and 


their  ranking  provides  a  basis  for  model  reduction.  In 
a  vibrations  application,  the  structural  modes  play  the 
part  of  "components,"  and  the  summed  contributions  of  all 
modes  make  up  the  value  of  the  total  system  response. 

Consider  the  matrix  second  order  system  in  Eq  (2)  and 
its  state  space  formulation.  In  the  case  of  small  damping, 
Eq  (7)  represents  a  set  of  dynamically  uncoupled  modal 
equations  in  the  open  loop  sense.  If  u  is  a  function  of 
state,  as  in  feedback  control,  these  equations  do  not  un¬ 
couple  (20) .  For  this  reason,  the  modal  cost  analysis 
techniques  are  split  into  open  and  closed  loop  forms.  This 
paper  will  deal  principally  with  the  open  loop  case. 

In  the  open  loop  case,  the  feedback  control  u  =  0  and 
we  are  interested  in  the  contribution  of  each  mode  to  the 
output  y  (open  loop  control  can  be  applied,  but  for  the 
purposes  here  the  modeling  of  open  loop  disturbance  inputs 
only  is  most  productive) .  A  cost  function  is  defined  as 

T  ^ 

V  =  lim  E(y  Qy}  =  EV.  (18 

t-°°  -  1  1 

that  is,  V  is  the  expected  mean  square  response  weighted 
by  the  matrix  Q,  and  is  made  up  of  the  sum  of  the  con¬ 
tributions  from  each  of  n  system  modes .  Note  that  y  here 
is  some  system  output  or  performance  criterion  of  interest; 
in  the  case  at  hand,  y  represents  the  deviation  from  a 
desired  pointing  constraint,  a  quantity  that  we  wish 


to  control. 


The  eigenvalues  of  the  damped,  unforced  system 


fi  + 

2C  •  co . 

n  +  j 

to2 

1  1 

p 

/ 

n  =  0  i=l . . . n 

are  given  by  the  roots  of  the  characteristic  equation 

=  0 

and  correspond  to  the  eigenvalues  of  the  matrix  A  in  Eq 
Their  form  is  given  by 
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(19) 


(20) 

(12)  . 


(21) 


where  the  approximation  holds  for  small  damping,  £^-*-0  i=l...n. 
In  addition,  the  eigenvalue  problem 


AE  =  EA 


for 

i=l , . . . 2n  all  (22) 

complex 
eigenvalues 


yields  the  matrix  E  whose  columns  are  the  so-called  "right" 
eigenvectors  of  A  corresponding  to  A^: 


E  = 


l  e 
•  “2  n 


(23) 


Here,  with  A  as  in  Eq  (16)  ,  E  takes  the  form 


The  "left"  eigenvectors  of  A  are  found  from 


FA  =  AF 

and  correspond  to  the  rows  of 


(25 


(26 


Now  define  an  "observability"  vector: 


=  Ce .  =  J 

[c  *' 

i  +  A  . 

\  c„$ 

1 

[-P-J 

k  i 

—V  — 

i=l , . . . 2n 
k=l , .  .  .n 
columns 


(27 


It  is  shown  in  (19,20)  that  the  individual  modal  cost 

for  the  open  loop  case  of  Eq  (12),  u  =  0,  can  be  expressed 

as 


V. 

l 


(28 


where  Q  is  the  weighting  matrix  from  the  cost  function, 

Eq  (18) ,  and  W  is  the  noise  strength  matrix  from  Eq  (3) . 
Note  that  if  the  disturbances  can  be  considered  independent 
and  uncorrelated,  W  is  a  diagonal  matrix  whose  elements  are 
the  mean  square  values  of  w(t) . 

To  show  this,  the  cost  V  may  be  written  in  terms  of 


4 


the  initial  condition  covariance  as 


V  =  E  /  (yTQy)  dt  =  tr(KP0) 


where 


K  =  /  e-  tCTgCe-t  dt 


PQ  =  E  x (0) ~x (0)  x (0)-x (0) 


Transforming  to  modal  coordinates  in  the  damped  system, 

Eq  (19) ,  the  elements  of  K  become,  in  terms  of  the  previous 
"observability"  vectors 

/oo 

e  (Xj.+Xj  )  t  dt  (32 

0 

=  -i/i^j  (2^Qgj 

T 

P0  transforms  as  PQ  =  EPE  ,  and,  from  Eq  (29), 

V  =  tr(KP0)  =  tr  (e“TKE-1EPET)  =  tr  (KP)  (33 


n  r 

=  Z  K.P. 
j=l-i-i 


i  =  ith 


column 


so  that 


v. =  kTp. 

i  -i-i 


(2iS2j)  (E^PQEj7)!  /(1,+Xj) 


where,  for  example,  E~  denotes  the  i  column  of  E- L . 

T  —  1  . 

Substituting  f^  *  E^  in  Eq  (34),  and  considering  the 
disturbance  statistics  rather  than  the  initial  condition 
statistics  (£«■*■  GWGT)  yields  Eq  (28). 


( 
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It  is  evident  from  Eq  (28)  that  in  the  general  case 
the  cost  due  to  an  individual  mode  is  the  result  of  the 


direct  excitation  of  that  mode  plus  the  coupling  effects 
of  all  other  modes.  However,  in  the  case  of  small  damping 
it  can  be  shown  (20)  that  the  term  resulting  from  i=j 
dominates  the  summation  to  the  extent  that  the  other  terms 
may  be  neglected.  This  is  equivalent  to  the  statement 
following  Eq  (7)  regarding  decoupling  of  a  lightly  damped 
system  in  modal  coordinates . 

Combining  the  cost  due  to  a  mode  and  its  conjugate 
in  Eq  (28)  leads  to  the  simplified  expression 


V.  = 
l 


(jQo.)f^(GWGT)fi 


Ai“i] 


(35) 


Denoting  the  ith  column  of  C  $  by  (C  <I>)  .  ,  the  terms  of  the 

-p-  -p-  i 

above  equation  become 

0 

oTQo.  =  (C  $)T  Q  (C  $)  .  +4)?  (C  <J>)T  Q  (C  <?)  .  +  0 A.  )  (36 

-1--1  -p-  1  -  -p-  1  1  -V-  1  -  -V-  1  jT  \ 


£i£WGTfi  =  (|TRWRT$)i  /  (4o)?) 


(37 


so  that  the  modal  costs  become 


V ,  =  < 


(c  $) :q(c  $) . 

-p-  1-  -p-  1 


“i(Svi,i2«Sv*)i 


T  T 
($  RWR1)  i 
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«'i“i 


(38 


ir  terms  of  explicitly  known  quantities  in  the  original 
system. 


In  the  study  here,  only  output  position  variables  were  nec¬ 


essary,  so  £v  =  0  and  the  expression  further  simplifies  to 

V,  =  (C  ®)T  Q  (C  $)  .  ($TRWRT$)  .  /(4c.w|)  (39) 

Thus,  based  on  the  system  model  and  its  computed  nat¬ 
ural  frequencies  and  mode  shapes ,  the  individual  modal 
contributions  to  an  "output"  can  be  assessed  and  ranked. 
Presumably,  modes  having  a  small  contribution  could  be 
truncated  from  a  control  model  whose  purpose  was  to  mini¬ 
mize  the  "output." 

A  similar  expression  can  be  developed  for  a  closed 
loop  system  with  u  /  0  in  Eq  (7) .  However,  this  case 
involves  feedback  gains  which  must  be  computed  first  for 
the  high  order  model;  closed  loop  modes  may  then  be  trun¬ 
cated.  The  computation  of  feedback  gains  for  a  high  order 
model  is  computationally  unwieldy,  to  say  the  least.  An 
alternative  approach  is:  first  apply  the  open  loop  method 

to  the  high  order  model,  compute  the  feedback  gains  for  a 
truncated  control  system,  and  then  apply  the  closed  loop 
version  of  the  model  cost  equation  for  further  model 
reduction.  No  closed  loop  applications  were  attempted 
here. 

Internal  Balancing 

A  third  model  reduction  technique  applicable  to  large 
systems  is  based  on  the  work  of  Moore  (14). 


Essentially,  Moore  proposes  that  a  system  representation 


which  has  been  transformed  to  an  "internally  balanced" 
coordinate  system  (Grammian  -  balanced)  may  be  state  - 
truncated  in  that  coordinate  system  to  produce  a  reduced 
order  model.  The  truncation  is  accomplished  based  on  the 
relative  controllability /observability  of  the  system 
states .  This  method  is  advantageous  in  that  it  incorpor¬ 
ates  information  about  sensor/actuator  location  that  is 
not  done  in  modal  cost  analysis;  in  addition,  any  chosen 
subset  of  the  transformed  system  model  is  guaranteed  to  be 
stable.  Unfortunately,  for  a  general  multi-input,  multi¬ 
output  system,  a  different  coordinate  basis  is  required 
for  each  input-output  pair.  However,  in  the  case  where 
damping  is  very  small,  it  can  be  shown  that  modal  coor¬ 
dinates  and  balanced  coordinates  are  approximately  equi¬ 
valent  (4,9).  Hence  the  technique's  applicability  to  the 
problem  at  hand. 

The  controllability  and  observability  Grammians 
(Gram  matrices)  derived  from  the  system  of  Eqs  (12),  (13), 
and  (14)  with  A  asymptotically  stable  and  time-invariant  are 


r- 

w.  =  /  e-tBBTe-”t 

J  Q 


dt 


(401 


f  00 

Wq  =  /  e-  tHTHe-t 

ft 


dt 


(41) 
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The  steady  state  Wc  and  WQ  satisfy  the  Lyapunov  matrix 
equations : 


AWq  +  WcAT  +  BBT  =  0 


ATWq  +  WqA  +  CTC  =  0 


(42) 


(43) 


The  eigenvalues  of  Wc  can  be . considered  to  be  a  qualitative 

measure  of  the  controllability  of  the  system  states  in  a 

given  coordinate  system  (17)  .  The  smaller  an  eigenvalue  of 

W  ,  the  more  difficult  it  is  to  affect  the  related  state 
-c 

by  an  application  of  control  energy;  i.e.,  it  is  less  con¬ 
trollable.  The  observability  Grammian  is  the  dual  of  Wc, 
and  hence  similar  properties  can  be  identified  in  Wq. 

It  is  desired  to  find  a  coordinate  system  wherein 
these  attributes  can  be  exploited  most  advantageously. 

Note  that  the  Grammians  are  not  invariant  under  coordinate 
transformations;  the  thrust  here  is  to  find  a  transformation 
that  will  leave  the  Grammians  equal  and  diagonal. 

The  system  of  Eqs  (12),  (13),  and  (14)  transforms  to 
this  new  coordinate  frame  through  the  transformation  T  as 


x  =  Tx 


A  =  T~ 1  AT 


C  =  CT 


B  =  T“ 1 B 


H  =  HT 


D  =  T“ 1 D 


(44) 


where 


bearing  in  mind  that  Uc  is  an  orthogonal  matrix  so  that 
- 1  T 

U  =  U  .  This  new  matrix  W*  decomposes  as 

No  “C  “  O 


w*  =  hth  =  u*E*2u*T 


-o 


-o-  -o 


where 


E *2  = 


\ 


a2 

i 


\ 


(49) 


(50) 


a  =  eigenvalues  of  W* 


U*  =  matrix  whose  columns  are  eigenvectors  of  W* 


It  is  shown  in  (9)  that  the  transformation 


T  =  U  Z^U*E*~is 
-  -c-c-o- 


(51) 


operates  on  W  and  WQ  such  that 


w  =  W  =  E* 
-c  -o 


(52) 


that  is,  the  controllability  and  observability  Grammians 
are  equal  and  diagonal  in  the  new  state  space.  As  well, 
the  diagonal  elements  of  E*2  are  the  eigenvalues  of  WqWc , 
and  the  elements  of 


E*  = 


\ 


a . 

i 


(53) 


are  singular  values  of  H  in  Eq  (4  8)  •  A  system  in  this 
form,  Wc  =  Wq  =  E*,  is  internally  balanced. 


i 
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Essentially,  then,  the  best  model  truncation  decisions 
based  on  controllability  and  observability  can  be  made 
by  rank  ordering  the  and  keeping  only  the  modes  corres¬ 
ponding  to  the  higher  values.  In  the  case  of  Eq  (45) , 
transformed  to  balanced  coordinates  and  partitioned  as 
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(54) 


with  >>  ,  then  input  u  affects  x^  less  than  it 

affects  x  ,  and  output  y  is  affected  by  less  than  by 

(4).  Note  that  additional  relationships  can  be  examined 
by  using  D  in  Eq  (40)  or  H  in  Eq  (41)  and  considering  the 
"disturbability "  and  "measureability "  Grammians. 

An  explicit  form  of  this  transformation  T  is  given  in 
(4) .  For  Eqs  (12)  and  (14)  with  w  =  0,  written  for  each 
mode  as 

(55) 


*i  =  dixi  +  Biu 
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■  I  Hi  x. 


for  kth  sensor 
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x :  i 


that  is, 
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=  kth  element  of  a  column 


tne  balancing  transformation  is 
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Applying  this  transformation  yields  a  balanced  subsystem 
with  balanced  Grammians 


n  - 


'n  ° 


In  the  case  of  small  damping,  c, ^  <<1,  then  Eqs  (60)  and  (61) 
reduce  to  =  Q 2\  ~  aif  an<*  t^ie  aPProxima'te  balanced 
system  is 


A.  x  .  +  B .  u 
-i-i  -i- 
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where 


x.  =  A.x. 


*i*i  Si  *  li  "Si  H1  =  HiTi 


-1  -T  T 

EJ  =  T.  AW„  T.  =  T7W  T. 
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The  major  result  here  (4)  is: 


1)  if  damping  is  small,  5^  <<  1,  then  the  ith  subsystem  of 
Eqs  (66)  -  (68)  is  approximately  internally  balanced 
with 


°i  = 


-1 


m  rn 

b .  b .  (n  .  h  .  + 
-1-1  -vi-vi 


(71) 


2)  if  all  modal  damping  is  small  and  the  modal  natural 
frequencies  are  sufficiently  distinct,  i.e.  if 


max(?i,Sj)max(<ui,a)j)/|ai  - 


CD 


<  < 


1,  i/j 


(72) 


then  the  entire  model  is  approximately  internally 
balanced  with 


Z*  = 


Z * 


i=l .  .  .n 


(73) 


Furthermore,  in  the  cases  of  position  or  velocity  sensors 

only,  that  is,  H  =  0  or  H  =0,  the  above  balanced  coor- 

p  v 

dinates  are  a  special  case  of  modal  coordinates.  In  other 
words,  no  transformation  is  needed  from  the  system  of 
Eqs  (12)  -  (14) ,  yet  the  above  results  still  apply.  As 
position  sensors  only  are  utilized  in  this  study,  this 
conclusion  is  valid. 

Several  cases  of  modal  ranking  may  now  be  investigated, 
depending  on  the  system  matrices  used  in  the  Grammian 
definitions,  Eqs  (40)  and  (41). 


Actuators  to  Output.  Being  consistent  with  the  system 
matrices  as  defined,  the  Grammians  of  interest  in  this  case 


The  balancing  approximation  yields  rankings  of  o__  ,  that 

BCi 

is,  the  degree  to  which  actuators  can  affect  the  ith  mode's 
contribution  to  the  output. 


Disturbances  to  Output.  One  may  consider  the  disturb¬ 
ance  propagation  through  the  system  to  the  output  by  replac¬ 
ing  B  with  G  in  Eq  (74).  The  resulting  a__  rankings  indicate 

GCi 

the  modes  that  are  most  necessary  to  control  in  order  that 
the  effects  of  disturbances  on  performance  is  minimized. 

Disturbances  to  Sensors.  Similarly,  replacing  C  by  H 
in  Eq  (75)  results  in  a  modal  ranking  of  the  relative  observ¬ 
ability  of  the  disturbance  inputs. 


Actuators  to  Sensors.  Using  the  B  and  H  matrices  in 
Eq  (74)  and  (75)  yields  a  ranking  of  the  ability  of  the 
actuator/sensor  configuration  to  sense  and  respond  to  a 
given  mode.  The  configuration  may  be  able  to  powerfully 
influence  a  mode,  but  if  this  mode  has  little  contribution 
to  the  performance,  this  ability  is  not  useful. 
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Utilizing  these  four  rankings,  one  may  analyze  the 


ft 


» 


ft 


modal  mechanism  of  the  disturbance  propagation  and  the 
corresponding  ability  of  the  given  actuator/sensor  con¬ 
figuration  to  sense  and  eliminate  these  disturbance 
effects.  This  provides  a  straightforward  method  for 
selective  truncation  of  the  high  order  mathematical  model 
to  the  low  order  control  design  model. 

A  relationship  exists  between  the  modal  costs  of 
Eq  (38)  and  the  internal  balancing  disturbance/output 
ranking  given  by 


'GCj.  = 


-1 


T  T  -2 

g .  g .  (c  .c  .  +  oj.  c  .  +  c  .) 
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(76) 


Assuming  disturbances  characterized  as  zero  mean  white 
noises  of  unit  strength ,  then 


V. 
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4  c  .  as  •  Q  ■  a_„ 
h  ri  GCi 


(77) 


Approximation  Error.  The  internal  balancing  approx¬ 
imation  is  valid  for  cases  that  satisfy  Eq  (72),  which  in¬ 
volves  modal  damping  ratios  and  natural  frequencies.  In 
this  study,  damping  is  assumed  uniform  for  all  modes  and 
equal  to  .005.  The  remaining  condition  requires  a  suffi¬ 
cient  separation  of  modal  natural  frequencies.  For  the 
case  of  ?  =  .005,  a  frequency  difference  of  5%  results  in 
an  approximation  error  of  less  than  5%  (4) .  Frequency  dif¬ 
ferences  of  less  than  this  result  in  significantly  higher 
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errors.  In  the  case  of  a  high  order  model,  modal  natural 
frequencies  may  in  fact  be  quite  close  together,  and  the 
approximation  error  may  be  unacceptable.  However,  the 
subset  of  modes  that  do  meet  the  error  criterion  still 
provides  useful  information  in  the  model  reduction  process; 
this  situation  will  be  discussed  in  connection  with  the 
CSDL  II  structure. 


IV.  CONTROL  ALGORITHM 


Once  an  open  loop  reduced  order  model  has  been  achieved, 
the  next  step  is  the  implementation  of  some  form  of  closed 
loop  (feedback)  control  technique.  In  this  study,  the  method 
utilized  is  that  developed  by  Calico  and  associated  students 
(1,3,8,23).  The  basic  methodology  is  presented  in  (8),  and 
will  be  summarized  here. 

State  Space  Formulation 

The  deterministic  feedback  control  system  is 


x  =  Ax  +  Bu  u  =  Gcx 

(78) 

y  =  Cx 

(79) 

z  =  Hx 

(80) 

This  system  is  in  the  form  of  the  mathematical  model  of 
Figure  1,  that  is,  it  is  a  finite  dimensional  representation 
of  the  infinite  dimensional  true  system.  The  modes  which 
have  been  eliminated  in  this  inherent  truncation  are  known 
as  unmodeled  modes;  they  are  typically  very  high  frequency 
and  do  not  significantly  impact  the  control  problem.  The 
modes  that  make  up  the  mathematical  model  of  Eqs  (78) -(80) 
are  the  modeled  modes. 

The  modeled  modes  are  identified  as  either  controlled, 
suppressed,  or  residual.  The  controlled  modes,  x_  are 

■"’G  f 

those  selected  by  the  model  reduction  techniques  that  must 
be  directly  controlled  in  order  for  the  system  to  meet  the 


performance  objectives.  The  remaining  modeled  but  uncon¬ 
trolled  modes  will  be  affected  by  the  application  of  con¬ 
trol  to  xc  even  though  they  themselves  are  not  directly 
controlled;  this  phenomenon  (see  Eq  (105) )  is  known  as 
control  spillover,  and  may  result  in  poor  performance  or 
even  closed  loop  instabilities  if  not  remedied.  The  un¬ 
controlled  modes  whose  reaction  to  this  spillover  must  be 

suppressed  are  designated  x_ ,  suppressed  modes.  The  re- 

s 

maining  modeled,  uncontrolled  modes  that  are  relatively 
insensitive  to  the  spillover  effect  are  termed  residual 
modes.  The  design  level  controller  is  unaware  of  the  exist- 
ance  of  these  residual  modes,  however  they  do  contribute 
to  the  performance  of  the  reduced  order  control  model;  hence 
their  inclusion  in  the  final  system. 

The  system  of  Eqs  (78) -(80)  may  be  separated  into  the 
states  corresponding  to  the  controlled,  suppressed,  and 
residual  modes  as 


x  =  A  x  +  Bu  u  =  Gx 

-c  -c-c  -c-  -  -c-c 


x  =  A  x  +  B  u 

-s  -s-s  -s- 


x  =  A  x  +  B  u 

-r  -r-r  -r- 


y  -  cx„  +  c  x„  +  c  x„ 

i  -c-c  -s-s  -r-r 


z  =  H  x  +  H  x  +  H  x 

-c-c  -s-s  -r-r 


The  unknown  quantity  here  is  the  feedback  gain  matrix 
G  •  This  is  found  by  the  application  of  steady  state  opti¬ 
mal  regulator  theory  (10) ,  where  a  quadratic  cost  function 
to  be  minimized  is 


J  = 
c 


T  m 

x  Q  X  +  u~R  u  dt 
-c-c-c  -  -c- 


Qc  s  ncxnc  state  weighting  matrix  (pos.  semi-definite) 
R  =  mxm  control  penalty  weighting  matrix  (pos.  def.) 


Here,  Rc  is  assumed  to  be  the  unit  diagonal.  The  optimum 

gain  G  to  minimize  the  cost  J  is  known  to  be 
c  c 


-1  T 

G  «  -R  lB'S 
— C  — c  -c— 


where  S  is  the  solution  to  the  steady  state  matrix  Riccati 
equation 


SA  +  ATS  -  SB  R_1BTS  +  Q  =  0 
— c  -c-  — c-c  -c-  sc 


However,  in  order  to  form  u  =  G  x  ,  the  state  x  is 

-  -c-c  -c 

required.  This  is  normally  not  known  explicitly,  but  only 
as  an  integral  part  of  the  measurement  vector  z  in  Eq  (85) • 

It  is  therefore  necessary  to  estimate  the  states  of  interest. 
In  this  deterministic  treatment,  a  Luenberger  observer  is 
introduced,  of  the  form 


x  =  Ax  +  Bu  +  K  ( z  -  z) 
-c  -c-c  -c-  -c  - 


'  -  *  ■  *  J  ' -A  ‘  *  '  *  *  •  ‘ -A  •  '  A  *  J*  *  J*  ’  A  ’  *  •  ■  •  •  • 


.  '-v  v-v 
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xc(0)  =  0 


z  = 

-  -c-c 


where  the  hat  notation  denotes  an  estimate.  Eq  (89)  results 
in  an  estimate  of  the  state  whose  error 


e  =  x  -  x 
-c  -c  -c 


approaches  zero  as  t-*-00  if  the  gain  Kc  is  chosen  so  that  the 
observer  is  asymptotically  stable.  By  combining  Eqs  (81), 
(85)  ,  (89) ,  and  (91) ,  the  estimator  can  be  written  in  terms 
of  the  error  states  as 


Sc  -  <AC  -  SCHC)SC  +  Kc(Hs£s  ♦  Hrxr) 


In  order  to  insure  the  desired  stability,  we  concentrate 
on  finding  the  observer  gain  matrix  Kc  such  that  the  eigen¬ 
values  of  A„  -  K  H  in 
-c  -c-c 


e  =  (A  -  KH)  e 
— c  -c  -c-c  -c 


retain  negative  real  parts  in  the  closed  loop  system.  To  be 
sure,  the  suppressed  and  residual  terms  in  Eq  (93)  still 
impact  the  observer  and  are  termed  observation  spillover 
(see  Eq  (105))  ;  the  method  of  suppressing  these  effects  will 
be  addressed  in  the  next  section. 


I 


In  order  to  determine  the  gain  K  ,  a  problem  analogous 
to  Eq  (94)  is  posed  as 


X  -  (*C  -  Me*  X 


•  „T 
w  =  A1  w 


where 


T 

g  =  K  w 
-  -c- 


Since  the  eigenvalues  of  a  matrix  equal  the  eigenvalues  of 
its  transpose,  a  gain  matrix  Kc  that  insures  stability  of 
Eq  (95)  will  also  satisfy  the  requirements  for  Eq  {94) . 

Eqs  (96)  and  (97)  respond  to  an  application  of  steady  state 
optimal  regulator  theory  as  before.  A  quadratic  observer 
cost  function  to  be  minimized  is 


f  r  T  t  l 

r  =  /  w  Q  w  +  g  R  g  dt 

o  J  I--O-  “°” J 


with  and  R  defined  in  a  manner  similar  to  those  in 
ao  -o 

Eq  (86)  .  These  weighting  matrices  should  be  chosen  so  that 
the  observer  error  in  Eq  (94)  exhibits  a  rapid  exponential 
decay.  The  optimal  gain  solution  of  Eq  (97)  is 


T  -1 

K  =  +R  H  P 
-c  -o  -c- 


where  P  is  the  solution  to  the  steady  state  matrix  Riccati 


equation 


PA^  +  A  P  -  PhJr^H  P  +0=0 
— c  -c-  — c-o  uc-  *0 


(100) 


1 


The  state  equations  (81)  ,  (82)  ,  and  (83;  may  now  be 
written  in  terms  of  the  gain  and  error  states  as 


x  =(A  +  BG)x  +  BGe 

—  C  — c  — *C“*C  “C  — c— G~C 


(101) 


5S  =  Ms  +  §sMc  +  §sMc 


(102) 


X  =  Ax  +  B  G  x  +  B  G  e 
— r  -r— r  -r-c-c  -r-c-c 


(103) 


Forming  the  augmented  state  vector  Z  as 


allows  the  expression  of  the  closed  loop  system  as 


(104) 


(*c  +  Sc5c> 


B  G 
-c-c 


(A-KH)  KH  KH 

-c  -c-c  -c-s  -c-r 


Z  (105) 


BG 

-s-c 


B  G 
-s-c 


B  G 

-r-c 


B  G 
-r-c 


4 


Note  that  the  application  of  control  gains  Gc  affects  not 

only  the  xc  modes  but  also  the  uncontrolled  modes  through 

the  control  spillover  terms  B  G  and  B  G  .  Likewise,  the 

— s-c  -r— c 

presence  of  the  suppressed  and  residual  modes  impacts  the 

observer  through  the  observation  spillover  terms  K  H  and 

KcHr.  The  effect  of  these  spillover  terms  may  be  small;  on 

the  other  hand,  as  control  and  observation  gains  G  and  K 

“C  "  c 

are  increased  in  an  effort  to  control  the  x  terms  more 

-c 

precisely,  the  spillover  terms  may  in  fact  cause  overall 
system  instabilities. 

Transformation  Technique 

In  light  of  the  above  mentioned  potential  for  insta¬ 
bilities,  a  method  is  desired  to  eliminate  or  suppress  the 
spillover  effects.  Eq  (105)  represents  the  full  mathematical 
system  model;  the  representation  of  only  those  modes  which 
require  active  intervention,  either  through  direct  control 
or  spillover  suppression,  is 


Z  = 


<5C  +  5cSc> 


B  G 
-s-c 


BG 

-c-c 


(A  -  K  H  )  K  H 

-c  -c-c  -c-s 


BG 

-s-c 


*s 


(106) 


where 


Z  = 


*c 


§c 


*s 


(107) 
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If  the  system  eigenvalues  of  Eq  (106)  can  be  guar¬ 
anteed  to  have  negative  real  parts,  then  this  system  will 
be  stable.  Now  the  diagonal  blocks  of  Eq  (106)  have  stable 

eigenvalues:  (A  +  B  G  )  and  (A  -  K  H  )  because  the  gains 

— c  — c~C 

have  been  formulated  so  that  this  is  so,  and  Ag  due  to  the 
natural  modal  damping.  Since  the  eigenvalues  of  a  block 
triangular  matrix  are  equal  to  the  eigenvalues  of  its  diag¬ 
onal  blocks,  a  block  triangularization  of  the  matrix  of 
Eq  (106)  would  produce  the  desired  stability.  Causing 
B  G  =  0  or  K  H  =  0  and  rearranging  accomplishes  this  tri- 
angularization.  Here,  a  transformation  is  sought  such  that 

KH  =  0. 

-c— s  - 

An  obvious  solution  is  K  =0;  however,  this  also 

-c  - 

renders  KcHc  *  0,  disrupting  the  observer  stability.  What 
is  needed  is  a  transformation  such  that 


KH  =  0 
-c-s  - 

(108) 

while 

K  H  t*  0 
-c-c  - 

(109) 

or,  in  the  form  convenient  to  the  following  development 


A  transformation  T  would  accomplish  this  if 


T 

H  T 

=  0 

(112) 

-s- 

T 

H  T 
-c— 

^  0 

(113) 

Assuming  such 

a  transformation  exists,  we  define 

a  vector 

2 

=  Tv 

(114) 

such  that  its 

substitution 

into  Eq  (96)  yields 

.T 

w  =  A  w  - 

T 

H  Tv 

(115) 

-c- 

-c — 

or 

,T 

w  =  A  w  - 

T 

H*  v 

(116) 

-c- 

-c  - 

with 

T  T 

H*  =  H^T 

(117) 

-c  -c 

and 

T 

v  =  K*  w 
-  -c  - 

(118) 

where  the  *  denotes  "transformed"  quantities.  Comparing 
Eqs  (115)  and  (118)  to  Eq  (96),  we  find 


T  T 

K  =  TK* 
-c  — c 


(119) 


However,  as  T  is  not  in  general  a  square  matrix,  the  new 
observer  gain  K*  is  found  as  before  from  the  solution  to  the 
steady  state  optimal  regulator  problem: 


K 


i< 


with  P  generated  from  the  steady  state  matrix  Riccati 
equation 


PAT  +  A  P  -  PH*TR*-1H*P  +  Q*  =  0 
— c  -c-  — c  -o  -c-  ao  - 


where 


R*  =  T  R  T 

-o  -  -o- 


(121) 


(122) 


Substituting  Eq  ( 119)  into  Eq  ( 110)  demonstrates  that  the 

gain  matrix  K*  in  fact  suppresses  the  observer  spillover 

terms  if  the  T  matrix  can  be  found  that  satisfies  Eqs  ( 112) 

and  ( 113) .  This  T  matrix  is  found  by  a  singular  value 

T 

decomposition  of  Hs . 

For  a  general  real  matrix  A+Rmxn,  the  following  is 
true  (22)  : 


T 

VAU  = 


diag  (a1  ,a2  , .  .  .ag) 


l  0 
0  0 


(123) 


1* 


U  2 


V 


singular  values  of  A  arranged  in  descending  order 
(square  roots  of  the  non-zero  eigenvalues  of  ATA) 

nxn  matrix  whose  columns  are  right  singular  vectors 
of  A  (eigenvectors  of  ATA) 

=  mxm  matrix  whose  columns  are  left  singular  vectors 
of  A  (eigenvectors  of  AAT) 

i  rank  of  A 


Thus,  the  number  of  non-zero  singular  values  equals  the 
rank  of  A.  Partioning  U  and  V  to  isolate  the  eigenvectors 
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i  — ' 


»  J 


.m*  .A  * 


i 


associated  with  the  non-zero  singular  values. 


u  =  u  u 

-q  .  -n-qj 

" 

V  =  VI  V  „ 
[-q  .  -m-qj 


(124) 


(125) 


From  Eq  (117) , 


a  =  v  nr 
-q — q 


(126) 


Since  U  is  an  orthogonal  matrix, 


UU  =  0 
-q-n-q 


(127) 


so  that 


AU  „  *  V  ZU'U  =  0 
— n-q  -q — q-n-q 


(128) 


In  the  problem  at  hand,  where  A  S  Hg ,  comparing  Eq  (104) 
and  Eq  (122),  the  required  transformation  is 


T  =  U 

-p-q 


(129) 


where 


=  matrix  whose  columns  are  the  eigenvectors  of  HgHs 
associated  with  its  zero  singular  values 

=  number  of  sensors 

2  rank  of  HT 
-s 


Assuming  H  to  be  a  full  rank  matrix  with  dimensions  pxnc 


in  the  case  of  position  sensors  only,  then  its  rank  will  be 
the  minimum  of  p  or  ns*  Eq  (129)  shows  that  in  order  that 
T  be  other  than  the  trivial  solution  T  =  0,  p  must  be 
greater  than  q;  that  is,  the  number  of  sensors  p  must  be 
greater  than  the  number  of  suppressed  modes  ng . 

A  parallel  development  to  eliminate  control  spillover 
terms  B  G  is  given  in  (1) ,  (8) . 
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V.  COMPUTER  IMPLEMENTATION 


Although  the  generation  of  a  working  computer  algor¬ 
ithm  was  an  extremely  time  consuming  portion  of  the  thesis 
process,  its  description  here  will  be  limited  to  brief 
comments.  A  source  code  listing  and  flow  diagram  may  be 
found  in  Appendix  B. 

The  computer  implementation  of  the  model  reduction 
techniques  was  fairly  straightforward.  The  code  was  de¬ 
signed  to  be  compatible  with  data  files  created  for  the 
controller  design  algorithm.  Input  consists  of  the  per¬ 
tinent  finite  element  analysis  results  for  the  model  of 
interest;  output  is  in  the  form  of  ranked  modes  from  the 
internal  balancing  scheme  and  from  the  modal  cost  analysis 
scheme.  In  addition  ,  a  check  is  made  of  the  appropriate¬ 
ness  of  the  approximate  internal  balancing  method  for  the 
case  at  hand. 

The  controller  design  and  spillover  suppression  algor¬ 
ithm  used  is  a  scaled  version  of  the  code  developed  by 
Aldridge  (1) .  In  its  original  form  the  code  was  designed 
for  the  investigation  of  multiple  controller  systems.  Here 
it  was  modified  for  use  with  one  controller  only,  as  the 
variable  of  interest  in  this  study  is  the  reduced  order 
model  and  not  the  control  scheme.  In  addition,  corrections 
were  made  to  the  time  response  portion  implemented  by  Baker 
(3)  so  that  true  system  output  versus  time  may  be  observed. 


Plotting  routines  and  the  use  of  procedure  files  allow 
the  output  of  eigenvalue  or  time  response  plots  via  the 
CDC  Cyber  NOS  interactive  system  with  only  one  user  in¬ 
put  system  command. 


t 


t 
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VI .  INVESTIGATION  AND  RESULTS 

The  CSDL  I  model  was  used  as  the  principle  test  case 
for  the  study.  Various  situations  were  investigated  involv¬ 
ing  the  different  model  reduction  schemes,  different  combin¬ 
ations  of  controlled  and  suppressed  modes,  and  different 
actuator/sensor  combinations .  Eigenvalue  plots  were  gener¬ 
ated  for  each  case;  in  addition,  plots  depicting  the  time 
response  of  node  one  under  the  action  of  the  various  feed¬ 
back  control  configurations  allowed  observation  of  the 
physical  reactions  to  each  arrangement.  These  time  response 
plots  track  the  radius  of  deflection  from  the  undeflected 
position  subject  to  the  initial  conditions  given  in  Appendix 
A.  In  all  cases  the  sampling  rate  is  20  Hz. 

Figure  4  shows  the  uncontrolled  response  and  the  effects 
of  controlling  the  first  four  modes  with  a  cost  function 
weighting  of  20.  It  is  desired  that  motion  due  to  the  initial 
conditions  be  effectively  damped  in  20  seconds;  this  time 
region  will  be  examined  in  subsequent  plots. 

Two  actuator/sensor  configurations  were  investigated: 

1)  six  co-located  actuator/sensor  pairs  on  the  "legs"  of  the 
structure,  and  2)  three  actuators  on  the  "arms"  of  the  tetra¬ 
hedron,  6  sensors  on  the  "legs"  (see  Figure  2).  Experimenta¬ 
tion  determined  that  a  control  weighting  of  20  made  the  re¬ 
sidual  modes  sufficiently  sensitive  to  spillover  to  be  ob¬ 
served  while  maintaining  reasonable  stability  in  the  unsup- 
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pressed  control  case.  Thus,  differences  due  to  the  model 
reduction  were  sought  in  the  unsuppressed  case;  often,  though, 
such  differences  were  not  apparent  until  suppression  was 
applied. 

Application  of  the  internal  balancing  and  modal  cost 
analysis  methods  yields  modal  rankings  for  the  CSDL  I  as 
given  in  Figure  5  and  6.  As  exactly  known  initial  condi¬ 
tions,  not  disturbances,  are  considered  for  the  CSDL  I, 
zero  mean  white  noise  disturbances  were  considered  to  be 
applied  at  the  actuator  locations  to  generate  the  modal . 
cost  rankings.  Weighting  the  output  variables  equally, 
the  modal  cost  rankings  are  essentially  the  same  as 
internal  balancing  except  for  the  inversion  of  modes  1 
and  2.  The  rankings  suggested  by  Figure  5  indicate,  for 
example,  that  although  mode  6  is  highly  controllable/observ¬ 
able  for  the  6/6  co-located  actuator/sensor  case,  its 
usefulness  in  achieving  the  desired  performance  objective 
is  minimal.  Thus  it  is  ranked  last  in  terms  of  modes  to 
be  controlled. 

Figure  7  is  included  as  an  aid  to  interpreting  the 
eigenvalue  plots.  In  all  cases,  weightings  for  the  con¬ 
trolled  mode  and  associated  observer  states  are  equal;  under 
this  condition  the  observer  eigenvalues  will  typically 
"shadow"  those  of  the  controlled  modes.  The  uncontrolled 
modes  remain  in  the  vicinity  of  their  naturally  damped  posi¬ 
tion,  disturbed  only  by  the  previously  discussed  spillover 


effects.  Note  that  mode  seven  is  driven  close  to  instability 
by  this  effect,  while  the  damping  of  mode  eight  is  actually 
improved.  Application  of  the  suppression  technique  of  Section 
IV  to  modes  6,  7,  and  8  removes  this  spillover;  however,  the 
tradeoff  required  is  evident  in  the  consequent  reduced  damp¬ 
ing  of  the  controlled  modes.  Note  that  mode  6  is  typically 
unaffected  by  spillover,  and,  although  easily  controlled, 
has  minimal  effect  on  the  line-of-sight  performance. 

Figures  8  thru  21  depict  various  controlled/suppressed 
mode  arrangements  chosen  according  to  frequency  truncation 
and  internal  balancing.  In  cases  exhibiting  instability 
caused  by  spillover,  extended  scale  time  responses  are 
included  to  show  the  long  term  system  behavior.  These 
instabilities  are  often  minimal  and,  in  concert  with  the  well 
damped  modes,  result  in  an  essentially  steady  state  motion 
as  shown  in  Figure  10.  More  severe  instabilities  are  obvious 
in  Figures  31  and  33.  The  degree  of  instability  for  a  given 
case  is  governed  by  the  cost  function  weightings  of  Eq  (86) 
and  (98).  Thus,  an  attempt  to  improve  the  output  response 
by  simply  adding  more  control  energy  to  the  controlled 
modes  may  have  the  opposite  effect. 

Unsuppressed  time  response  differences  brought  about 
solely  due  to  the  choice  of  model  reduction  schemes  are 
not  significant  until  the  case  of  5  controlled  modes  in 
Figures  18  and  20.  Here,  the  best  performance  before  appli¬ 
cation  of  suppression  was  achieved  by  controlling  5  modes 
according  to  the  balancing  selection  process  as  in  Figure  20. 


Suppression  actually  degraded  the  performance  somewhat,  as 
the  damping  of  mode  3  had  been  enhanced  by  the  spillover; 
elimination  of  this  spillover  reduced  the  damping  and  di¬ 
verted  control  energy  from  the  controlled  modes  to  accom¬ 
plish  this. 

A  conclusion  that  might  be  reached  is  that  one  could 
examine  the  controlled,  unsuppressed  eigenvalues  and  suppress 
only  those  modes  whose  damping  was  degraded  by  spillover. 
However,  a  test  case  of  suppressing  only  mode  9  in  Figure  20 
resulted  in  poorer  performance,  as  mode  3  was  driven  unstable 
in  the  process.  A  predictable  relationship  between  spill¬ 
over  and  mode  stability  is  elusive. 

In  addition,  there  are  situations  when  the  balancing 
algorithm  chooses  a  mode  to  be  suppressed  that  has  been 
driven  unstable  by  spillover,  whereas  simple  frequency 
truncation  does  not.  This  is  true  when  controlling  only 
the  two  most  critical  modes  as  shown  in  Figures  8  thru  13; 
time  responses  after  suppression  reflect  this  difference. 

A  second  actuator/sensor  configuration  was  examined 
consisting  of  three  actuators  located  on  the  tetrahedron 
arms  and  six  sensors  on  the  legs.  Such  a  setup  represents 
a  possible  controller  designed  to  reduce  hardware  complexity 
and  computation  requirements.  Corresponding  internal  bal¬ 
ancing  rankings  and  modal  costs  are  given  in  Figures  22  and 
23.  Figures  29  thru  33  illustrate  some  representative 
controller  behaviors  for  frequencv  truncated  and  internally 
balanced  reduced  models.  An  interesting  effect  is  demon- 


strated  in  the  unsuppressed  time  responses  of  Figures  25 
and  28.  Although  both  are  unstable,  the  balancing  model 
is  much  more  so,  illustrating  the  fact  that  the  choice  of 
a  reduced  order  model  by  the  one  criterion  of  actuator/out- 
put  modal  ranking  does  not  predict  what  the  closed  loop  be¬ 
havior  will  be.  Here,  the  balancing  algorithm  does  not  an¬ 
ticipate  the  instability  introduced  by  the  spillover  effect 
on  mode  4.  When  suppression  is  incorporated,  however,  the 
balanced  model  outperforms  that  generated  by  frequency  trunc¬ 
ation. 

Internal  balancing  modal  rankings  for  the  first  thirty- 
three  modes  of  the  CSDL  II  model  are  given  in  Figures  34  and 
35.  Assuming  unit  intensity  white  noise  disturbances  at  the 
actuator  locations,  modal  costs  in  Figure  36  yield  identical 
rankings.  The  first  three  modes  are  attitude  rigid  body 
modes;  they  are  immediately  included  in  any  controller  design. 
The  modal  observability /controllability  given  by  the  actuator/ 
sensor  ranking  matches  reasonably  well  with  the  modes  most 
significant  to  the  output;  it  can  be  expected  from  this  that 
the  actuator/sensor  configuration  would  be  able  to  damp  those 
modes  most  critical  to  the  performance.  Time  responses  for 
the  CSDL  II  model  were  not  generated  in  this  study. 

As  discussed  previously,  for  higher  order  models  the 
distinctness  of  the  natural  frequencies  may  present  a  problem. 
Modes  for  which  the  approximate  internal  balancing  algorithm 
is  seriously  in  error  are  given  in  Table  1.  For  the  most 
significant  modes  as  ranked  in  Figures  34  and  35,  only  mode  6 


does  not  meet  the  distinctness  requirement  for  a  valid  ap¬ 
proximate  balanced  model.  Thus,  deletion  of  all  the  modes 
in  Table  1  still  leaves  a  usable  ranking  of  modes  most 
critical  to  the  control  of  the  CSDL  II  line-of-sight  point¬ 
ing  error. 


RELATIVE  COST 


MODAL  COST 


■  J 


FIGURE  6 .  Modal  Cost  Rankings,  CSDL  I,  6/6  act/sen 
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FIGURE  16.  Internal  Balancing,  CSDL  I,  6/6  act/sen 
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FIGURE  19 .  Time  Response  -  Frequency  Truncation,  CSDL  I 
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FIGURE  25.  Time  Response  -  Frequency  Truncation,  CSDL  I 
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FIGURE  27.  Internal  Balancing,  CSDL  I,  3/6  act  sen 
4  controlled  -  2  suppressed  modes 
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FIGURE  30.  Time  Response  -  Frequency  Truncation,  CSDL  I 

4  controlled  -  4  suppressed  modes;  3/6  act/sen 
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FIGURE  35.  Internal  Balancing  Rankings,  onH 
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FIGURE  36.  Modal  Cost  Rankings,  CSDL  II,  21/21  act/sen 
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VII.  CONCLUSIONS 


The  approximate  internal  balancing  method  offers  a 
conceptually  sound  method  for  the  choice  of  a  reduced  order 
control  model.  This  method  utilizes  information  about  the 
sensor/actuator  configuration  in  the  mode  selection  pro¬ 
cess;  modal  cost  analysis  does  not.  In  the  CSDL  I  test 
cases  it  was  not  possible  to  demonstrate  a  clear  super¬ 
iority  of  one  method,  since  the  modes  most  critical  to  the 
performance  objective  (1  and  2)  were  contained  in  all 
reduced  models.  However,  the  improvement  in  time  response 
seen  in  Figure  18  would  seem  to  suggest  the  use  of  the 
internally  balanced  model  for  that  case. 

The  3  actuator/6  sensor  test  cases  indicate  that  com¬ 
parable  performance  may  be  achieved  by  a  reduction  in  the 
number  of  actuators  on  CSDL  I.  However,  spillover  suppres¬ 
sion  must  be  applied  to  eliminate  instabilities. 

Application  of  the  methods  to  the  CSDL  II  model  gen¬ 
erated  modal  rankings  that  are  significantly  different  from 
a  simple  frequency  truncation.  The  internal  balancing  method 
reveals  modal  contributions  to  the  system  performance  that 
are  not  otherwise  obvious,  and  indicates  the  ability  of 
the  sensor/actuator  configuration  to  influence  these  modes. 
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VIII.  RECOMMENDATIONS 


The  following  recommendations  are  suggested  for  follow- 
on  research: 

1)  .  The  effect  of  the  reduced  order  model  on  the  sta¬ 

bility  and  time  response  of  the  CSDL  II  should  be 
investigated. 

2)  .  It  was  assumed  here  that  ideal  sensor/actuators 

existed.  The  validity  of  this  assumption  could  be 
investigated  by  modeling  the  actuator/sensors  as 
noisy  devices  and  approaching  the  control  problem 
from  a  stochastic  viewpoint. 

3)  .  Durability  will  be  of  concern  in  an  actual  system. 

It  would  be  useful  to  examine  the  control  system 
robustness  in  the  face  of  sensor/actuator  failures. 

4)  .  As  selection  of  a  reduced  order  model  and  controller 

design  is  so  dependent  on  the  accuracy  of  the  struc¬ 
tural  model,  it  would  be  productive  to  investigate 
an  adaptive  control  law  that  would  account  for 
modeling  inaccuracies  and  changes  in  the  structure 
from  earth  to  space  environments.  The  choice  of 
modes  to  be  controlled  and  suppressed  could  possibly 


BIBLIOGRAPHY 


1.  Aldridge,  E.S.  Decentralized  Control  of  a  Large 
Space  Structure  as  Applied  to  the  CSDL  2  Model. 

MS  Thesis.  Wright-Patterson  AFB ,  Ohio:  School 

of  Engineering,  Air  Force  Institute  of  Technology, 

December  1982. 

2.  Balas,  M.J.  "Trends  in  Large  Space  Structure  Control 

Theory:  Fondest  Hopes,  Wildest  Dreams,"  IEEE  Trans. 

Auto .  Control ,  Vol .  AC-27,  No. 3,  1982. 

3.  Baker,  R.G.  Investigation  of  Multiple  Controller 
System  Efficiency  as  Applied  to  the  CSDL  I  Model. 

MS  Thesis.  Wright-Patterson  AFB,  Ohio:  School  of 
Engineering,  Air  Force  Institute  of  Technology, 

March  1984. 

4.  Gregory,  C.Z.  "Reduction  of  Large  Flexible  Spacecraft 
Models  Using  Internal  Balancing  Theory,"  AIAA  Guidance 
and  Control  Conference,  Gatlinberg,  Tenn. ,  1983. 

5.  Hablani,  H.B.  "Synthesis  of  a  Modern/Classical 
Controller  for  Large  Space  Structures , "  AIAA  Guidance 
and  Control  Conference,  San  Diego,  Ca.,  1982. 

6.  Hanks,  B.R.  and  Skelton,  R.E.  "Designing  Structures  for 
Reduced  Response  by  Modern  Control  Theory,"  Proceedings  of 
AIAA/ASME/ASCE/AHS  24th  Structures,  Structural  Dynamics, 
and  Materials  Conference,  Lake  Tahoe,  Nevada,  1983. 

7.  Hughes,  P.C.  and  Skelton,  R.E.  "Controllability  and 
Observability  for  Flexible  Spacecraft,"  AIAA  Journal  of 
Guidance  and  Control.  Vol.  3,  No. 5,  1980. 

8.  Janiszewski,  A.M.  Modern  Optimal  Control  Methods 
Applied  in  Active  Control  of  a  Tetrahedron.  MS  Thesis. 
Wright  Patterson  AFB,  Ohio:  School  of  Engineering,  Air 
Force  Institute  of  Technology,  December  1980. 

9.  Jonckheere,  E.A.  and  Silverman,  L.E.  "Singular  Value 
Analysis  of  Deformable  Systems , "  Proceedings  of  IEEE 
Conference  on  Decision  and  Control,  San  Diego,  Ca., 
December  1981,  pp.  660-668. 

10.  Kwakernaak,  H.,  and  Sivan,  R.  Linear  Optimal  Control 
Systems .  New  York:  John  Wiley  &  Sons  Inc.,  1972. 

11.  Laub,  A. J.  "Computation  of  'Balancing'  Transformations," 
Proceedings  of  Joint  Automatic  Control  Conference,  San 
Francisco,  Ca.,  August  1980. 


4 


87 


r 


% 


« 


12.  Meirovitch,  L.  Elements  of  Vibration  Analysis.  New 
York:  McGraw-Hill  Book  Company,  1975. 

13.  Moore,  B.C.  "Singular  Value  Analysis  of  Linear  Systems," 
Proceedings  of  IEEE  Decision  and  Control  Conference,  1978. 

14.  Moore,  B.C.  "Principal  Component  Analysis  in  Linear 

Systems:  Controllability,  Observability,  and  Model 

Reduction,"  IEEE  Trans.  Auto.  Control,  Vol.  AC-26, 

No . 1 ,  February  1981. 

15.  Pernebo,  L.  and  Silverman,  L.M.  "Model  Reduction  via 
Balanced  State  Space  Representations,"  IEEE  Trans .  Auto . 
Control ,  Vol.  AC-27,  No. 2,  April  1982. 

16.  RADC-TR- 82-21.  ACOSS  FIVE  (Active  Control  of  Space 
Structures)  Phase  IA.  Study  conducted  by  Lockheed 
Missiles  and  Space  Company,  Inc.  for  the  Rome  Air  De¬ 
velopment  Center,  March  1982  (ADA116655) . 

17.  Reid,  J.G.  Linear  System  Fundamentals .  New  York: 
McGraw-Hill  Book  Company,  1983. 

18.  Skelton,  R.E.  "Cost  Decomposition  of  Linear  Systems 
With  Application  to  Model  Reduction,"  International 
Journal  of  Control,  Vol.  32,  No. 6,  1980. 

19.  Skelton,  R.E.  and  Gregory,  C.Z.  "Measurement  Feedback  and 
Model  Reduction  by  Modal  Cost  Analysis,"  Proceedings , 

Joint  Automatic  Control  Conference,  Denver,  Colorado,  1979. 


20.  Skelton,  R.E.  and  Hughes,  P.C.  "Modal  Cost  Analysis  for 
Linear  Matrix-Second-Ordev  Systems , "  Journal  of  Dynamic 
Systems,  Measurement,  and  Control,  Vol.  102/151,  September 
1980. 

21.  Skelton,  R.E.  and  Hughes,  P.C.  and  Hablani,  H.B.  "Order 
Reduction  for  Models  of  Space  Structures  Using  Modal  Cost 
Analysis,"  Journal  Guidance,  Control,  and  Dynamics,  Vol. 
5,  No. 4,  July-August  1982. 

22.  Stewart,  G.W.  Introduction  to  Matrix  Computations. 

New  York:  Academic  Press,  Inc.,  1973. 

23.  Thyfault,  D.V.  Decentralized  Control  of  a  Large  Space 
Structure  Using  Direct  Output  Feedback.  MS  Thesis. 
Wright-Patterson  AFB ,  Ohio":  School  of  Engineering, 

Air  Force  Institute  of  Technology,  December  1983. 


4 


4 


88 


AD-A151  784  MODEL  REDUCTION  TECHNIOUES  APPLIED  TO  THE  CONTROL  OF  2/2 

LARGE  SPACE  STRUCTURES(U)  AIR  FORCE  INST  OF  TECH 
HRIGHT-PATTERSON  AFB  OH  SCHOOL  OF  ENGI.  .  D  U  VRRHOLA 
UNCLASSIFIED  DEC  84  AFIT/GA/AA/84D-11  F/G  14/2  NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  0f  STANDARDS- 1963- A 


< 


Appendix  A 


1 

l‘ 

) 

► 

CSDL  I  Node  Coordinates 

■  ■■  m  m  ■  w— ■  +  ,  »  ~  w  «.'  j  yj 

i 

• 

Node 

X  Y 

2 

§ 

1 

0.0  0.0 

10.165 

2 

-5.0  -2.887 

2.0 

3 

5.0  -2.887 

2.0 

4 

0.0  5.7735 

2.0 

‘  . 

•  ,  ,  •*» 

5 

-6.0  -1.1547 

0.0 

J  -  *** 

6 

-4.0  -4.6188 

0.0 

» 

7 

4.0  -4.6188 

0.0 

8 

6.0  -1.1547 

0.0 

9 

2.0  5.7735 

0.0 

■- 

10 

-2.0  5.7735 

0.0 

• 

• 

c 

CSDL  I  Initial  Conditions 

i 

Mode 

Displacement  (n) 

Velocity  (n) 

1 

-.001 

-.003 

? 

2 

.006 

.01 

3 

.001 

.03 

• 

4 

-.009 

-.02 

5 

.008 

.02 

6 

-.001 

-.02 

,  r 

7- 

-.002 

-.003 

i 

8 

.002 

.004 

9 

.0 

.0 

1  ( 

10 

.0 

.0 

i 

11 

.0 

.0 

12 

.0 

.0 

c 

1 

90 

1  ( 

1 

1 

SSBBfl 

CSDL  I  D  Matrix  (Actuator  Map) 
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NASTRAN  Eigenvalue  Analysis  Results  for  CSDL  I  (Nominal) 


Mode 
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Mass 
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Stiffness 

rad 
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7.29E+01 

8 . 539E+00 

7 . 29E+01 

10 

1.0E+00 

8.56E+01 

9 . 250E+00 

8.56E+01 

11 

1.0E+00 

1.05E+02 

1.029E+01 

1.05E+02 

12 

1.0E+00 

1.67E+02 

1 . 291E+01 

1.67E+02 

1 

ft 

• 

CSDL 

II  Node 

Coordinates 

(Meters) 

• 

• 

Node 

X 

Y 

Z 

Node 

X 

Y 

Z 

1 

-7.0 

0.0 

0.0 

37 

-4.0 

-3.0 

24.0 

• 

2 

-4.0 

5.0 

0.0 

38 

4.0 

3.0 

24.0 

ft 

3 

-4.0 

-5.0 

0.0 

39 

4.0 

-3.0 

24.0 

4 

0.0 

5.0 

0.0 

40 

0.0 

2.5 

2.0 

5 

4.0 

5.0 

0.0 

42 

0.0 

5.0 

-3.0 

# 

6 

4.0 

-5.0 

0.0 

43 

-2.0 

0.0 

-1.3 

ft 

7 

7.0 

0.0 

0.0 

44 

0.0 

-1.7 

-1.3 

8 

-7.0 

0.0 

2.0 

45 

2.0 

0.0 

-1.3 

9 

-4.0 

5.0 

2.0 

46 

-4.0 

-5.0 

-0.3 

« 

10 

-4.0 

-5.0 

2.0 

47 

4.0 

-5.0 

-0.3 

ft 

11 

4.0 

5.0 

2.0 

48 

-26.0 

0.0 

-1.3 

12 

4.0 

-5.0 

2.0 

49 

21.0 

0.0 

-1.3 

13 

7.0 

0.0 

2.0 

50 

16.0 

0.0 

-1.3 

ft 

14 

-6.0 

0.0 

12.0 

51 

-11.0 

0.0 

-1.3 

ft 

15 

-4.0 

4.0 

12.0 

52 

-6.0 

0.0 

-1.3 

16 

-4.0 

-4.0 

12.0 

53 

6.0 

0.0 

-1.3 

17 

4.0 

4.0 

12.0 

54 

11.0 

0.0 

-1.3 

’■■■■;  •' ' 

• 

18 

4.0 

-4.0 

12.0 

55 

16.0 

0.0 

-1.3 

ft  ' 

19 

6.0 

0.0 

12.0 

56 

21.0 

0.0 

-1.3 

26 

-5.0 

0.0 

22.0 

57 

26.0 

0.0 

-1.3 

27 

-4.0 

3.0 

22.0 

100 

0.0 

0.0 

0.0 

4 

28 

-4.0 

-3.0 

22.0 

910 

-4.0 

-2.5 

2.0 

• 

29 

4.0 

3.0 

22.0 

1001 

0.0 

-6.5 

22.0 

*  .* 

30 

4.0 

-3.0 

22.0 

1002 

0.0 

0.0 

2.0 

31 

5.0 

0.0 

22.0 

1003 

0.0 

6.5 

22.0 

• 

32 

-4.0 

10.0 

22.0 

1004 

0.0 

4.0 

2.0 

• . 

33 

4.0 

10.0 

22.0 

1112 

4.0 

-2.5 

2.0 

•  .  ■** . 

34 

-4.0 

-10.0 

22.0 

2830 

0.0 

-3.0 

22.0 

35 

4.0 

-10.0 

22.0 

3233 

0.0 

10.0 

22.0 

j 

« 

•  J 

35 

-4.0 

3.0 

24.0 

1 

« 
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Node 

CSDL  II  Lumped  Mass 

Distribution 

Mass  (k< 

Mass  (kg) 

Node 

9 

67.4 

44 

3500 

10 

67.4 

48 

81.9 

11 

67.4 

50 

1,638 

12 

67.4 

52 

73.8 

27 

69.5 

53 

73.8 

28 

6.74 

55 

163.8 

29 

69.5 

57 

81.9 

30 

6.74 

1001 

1000 

32 

6.74 

1002 

800 

33 

6.74 

1003 

1200 

34 

69.5 

1004 

600 

35 

69.5 

CSDL  II  Sensor/Actuator  Locations  and  Orientations 


i 


NASTRAN 

Eigenvalue 

Analysis  Results  for  CSDL 

II  (Nominal) 

• 

Mode 

Generalized 

Mass 

Generalized 

Stiffness 

rad 
w  sec 

a  ZSl-l 

sec 

1-6 

1.00E+00 

0.0 

0.0 

0.0 

7 

1.00E+00 

5 . 128E-01 

7.161E-01 

5.128E-01 

• 

8 

1.00E+00 

8 . 521E-01 

9 . 231E-01 

8.521E-01 

9 

1 . 00E+00 

8 . 835E-01 

9.399E-01 

8 . 835E-01 

10 

1.00E+00 

1.212E+00 

1 . 101E+00 

1.212E+00 

11 

1.00E+00 

8 . 819E+00 

2 . 862E+00 

8 . 189E+00 

• 

12 

1.00E+00 

1.266E+01 

3.502E+00 

1 . 226E+01 

13 

1.00E+00 

1 . 403E+01 

3 . 746E+00 

1 . 403E+01 

14 

1.00E+00 

1 . 492E+01 

3 . 863E+00 

1 . 492E+01 

15 

1.00E+00 

1 . 599E+01 

3.998E+00 

1 . 599E+01 

< 

16 

1.00E+00 

1 . 625E+01 

4 . 032E+00 

1.625E+01 

17 

1.00E+00 

2 . 623E+01 

5.122E+00 

2.623E+01 

18 

1.00E+00 

2 . 630E+01 

5.128E+00 

2 . 630E+01 

19 

1.00E+00 

2 . 677E+01 

5 . 17  4E+00 

2.677E+01 

• 

20 

1.00E+00 

3 . 310E+01 

5.753E+00 

3 . 310E+01 

21 

1.00E+00 

3 . 730E+01 

6 . 197E+00 

3.730E+01 

22 

1.00E+00 

5 . 301E+01 

7.281E+00 

5.301E+01 

23 

1.00E+00 

9 . 489E+01 

9 . 746E+00 

9.498E+01 

• 

24 

1.00E+00 

1.241E+02 

1.114E+01 

1.241E+02 

25 

1.00E+00 

1 . 999E+02 

1.414E+01 

1 .999E+02 

26 

1.00E+00 

2 .  001E+02 

1.416E+01 

2 . 001E+02 

■  „ 

27 

1.00E+00 

4 . 654E+02 

2.157E+01 

4 -654E+02 

i 

28 

1.00E+00 

4 . 705E+02 

2.169E+01 

4 . 705E+02 

29 

1.00E+00 

6 . 182E+02 

2 . 468E+01 

6 . 182E+02 

30 

1.00E+00 

6 . 276E+02 

2.505E+01 

6.275E+02 

31 

1.00E+00 

6 . 481E+02 

2.546E+01 

6.481E+02 

i 

32 

1.00E+00 

7 . 428E+02 

2 . 725E+01 

7 . 428E+02 

33 

1.00E+00 

1.700E+03 

4.123E+01 

1.700E+03 

34 

1.00E+00 

2 . 568E+03 

5 . 067E+01 

2 . 568E+03 

35 

1.00E+00 

2 . 821E+03 

5.311E+01 

2.821E+03 

!  C 

36 

1.00E+00 

3 . 095E+03 

5 . 563E+01 

3 . 095E+03 

y 


xlos 

YLOS 

ZLOS 


modal  coordinates* 


y  = 


n  = 


- . 6234E-13 

. 7492E-02 

. 1045E-08 

- . 7284E-05 

. 1051E-05 

. 1453E-12 

- .  5364E-03 

. 8123E-07 

. 1103E-13 

. 1777E-02 

- . 1218E-0  5 

- . 3412E-05 

- . 1880E-05 

- . 4005E-04 

. 3185E-05 

- .  7555E-08 

- . 1351E-02 

- . 8232E-06 

- .  9693E-03 

- . 1106E-04 

. 1901E-05 

-.5472E-04 

- . 8429E-04 

. 8500E-04 

. 6169E-02 

. 2150E-04 

. 1431E-04 

. 7284E-07 

- . 7115E-02 

- . 8065E-04 

- . 4461E-02 

- . 507  IE-04 

. 3427E-05 

. 2215E-04 

- . 1870E-02 

- . 6674E-04 

- . 6301E-03 

. 6725E-04 

- . 1256E-04 

- . 1051E-04 

- . 2819E-03 

. 76 12E-05 

. 4947E-05 

-.8486E-03 

- . 2771E-04 

- . 7761E-15 

- . 2288E-12 

- . 1626E-15 

- . 2869E-04 

- . 2736E-03 

. 3586E-04 

- . 3600E-02 

- . 8429E-04 

. 1385E-04 

- . 2285E-01 

- . 1862E-05 

- . 4  502E-04 

. 1621E+00 

. 1772E-03 

. 2510E-03 

. 6871E-02 

. 3786E-03 

- . 1182E-03 

- . 2481E-14 

- . 8148E-12 

- . 3790E-16 

. 4950E-06 

. 2705E-04 

- . 3122E-05 

. 2719E-13 

- . 3934E-11 

. 2880E-17 

- . 1781E-05 

- . 4463E-04 

. 5171E-05 

- . 1051E-05 

- . 5921E-04 

- . 1644E-04 

. 3466E-05 

. 1669E-04 

- . 2214E-04 

- . 3195E-03 

. 6  525E-05 

. 2638E-06 

- . 1086E-02 

. 1486E-05 

- . 2936E-06 

- . 4526E-03 

1765E+00 

- .  6466E-01 

. 1047E+00 

. 9565E-02 

- .  6518E-02 

. 1313E-01 

. 1605E-01 

.  5918E-01 

- . 7264E-02 

- . 1778E-02 

. 1289E-01 

*  here,  3  attitude  rigid  body  modes  and 
30  flexible  modes 


Appendix  B 
Source  Code  Listing 
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